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Abstract. 

We suggest that superscaling analyses of few-GeV inclusive electron scat- 
tering from nuclei, both in the quasielastic peak and in the region where the 
A-excitation dominates, allow one to make reliable predictions for charge- 
changing neutrino reactions at energies of a few GeV, relevant for neutrino 
oscillation experiments. 



1 Introduction 

Neutrino scattering has been the object of several recent investigations in con- 
nection with ongoing and planned experiments exploring neutrino properties, 
which use nuclei as targets. Moreover, the availability of new neutrino beams 
and nuclear targets opens the possibility of extracting information on the nu- 
cleon's properties, in particular on its strange form factors. For both purposes 
good control of the nuclear effects is essential for a reliable interpretation of the 
data. 

Present analyses of neutrino oscillation experiments make use of the Rel- 
ativistic Fermi Gas model, which, although reproducing the gross features of 
inclusive electron scattering responses, does not account for their detailed struc- 
ture. As will be shown in the following, relativitistic effects play an important 
role at the kinematical conditions relevant for modern experiments, namely neu- 
trino energies of a few GeV: hence the necessity of more realistic nuclear mod- 
eling, which includes relativistic dynamics. 
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While a number of relativistic calculations for neutrino scattering have been 
performed recently [1-7], a different approach can also be taken [8,9]: since any 
reliable model for neutrino-nucleus cross sections should first be tested against 
electron scattering experimental data, one can try to extract the neutrino cross 
sections from the experimental (e, e') cross sections, where a large amount of 
data is available. 

It has been shown [8,9] that, in appropriate kinematical conditions, this start - 
egy can be pursued as a consequence of the superscaling properties of the elec- 
tron scattering data. We shall illustrate this procedure after a brief review of the 
formalism. 

2 Electron and neutrino scattering: formalism and Relativistic 
Fermi Gas 

Electron and neutrino scattering off complex nuclei are closely related processes 
and as such they can be treated within the same formalism. 

Apart from obvious differences in the kinematics, due to the leptonic masses 
involved in the two processes, the basic difference between e-A and v-A scat- 
tering is the nature of the exchanged vector boson, a virtual photon, probing the 
electromagnetic nuclear current, in the former case and a W ± or a Z°, probing 
the weak current, in the latter. 

As a consequence, the corresponding cross sections involve different nuclear 
response functions and can be in general written in a Rosenbluth-like form as: 
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for inclusive electromagnetic (in the extreme relativistic limit m e — > 0), weak 
charged-current (CC, where the outgoing lepton I is detected) and weak neutral 
current (NC, where the nucleon N is knocked out) reactions, respectively. In the 
above (JMott the Mott cross section and a^ CC ' NC ^ the analogous quantities for 
the weak processes, Vi and Vi are factors depending upon the lepton kinematics 
(see Refs. [8, 9] for their explicit expressions) and the ± sign refers to neutrino 
or antineutrino scattering. 

Note that whereas the (e, e') cross section depends on two electromagnetic 
respone functions, Rl.t, three weak responses, Rl,t,t', enter the CC neutrino 
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scattering due to the presence of the axial weak current. In the NC reaction 
three more responses, Rtt,tl,tl' are involved, similarly to what happens in 
the semi-inclusive (e, e'N) process. 

Let us first focus on the quasielastic peak, where the dominant process is the 
knockout of a single nucleon. This domain can be treated in first approximation 
within the relativistic Fermi gas (RFG) model, where the nucleus is described 
as a collection of non-interacting Dirac particles. In this case all the response 
functions in ( I1I2I3> turn out to be given by 

(i? 4 )RFG = ^LRi*^ /rF G (V0> (4) 

qkp 

where TV is the appropriate nucleon number, tojv the nucleon mass, q the three- 
momentum transfer and kp the Fermi momentum. Moreover Ri S n ' are the 
corresponding single nucleon responses and 

/rfgW0 = |(i-V 2 Mi-V> 2 ) (5) 

is the RFG superscaling function, which only depends upon a specific combina- 
tion of the variable q, to (energy transfer) and kp, namely the scaling variable 

^ = ±^/J, T = (ly/i + l/T-w/2-m N ). (6) 

In Eq. Tp is the Fermi kinetic energy and To the minimum energy a nucleon 
in the Fermi sphere must have in order to participate in the reaction at fixed q and 
oj, t = {q 2 —uj 2 )/ (Am 2 N ) being the dimensionless four- momentum transfer. The 
± sign in refers to values of uj smaller (— ) or larger (+) than the quasielastic 
peak position uj = \J q 2 + m' 2 N — m^, corresponding to i)j = 0, and the response 
region is limited by the ^-function in (|5} to the range — 1 < < 1. 

We stress that all the response functions, and hence the cross section, can 
be factorized in a "single-nucleon" function times the superscaling function /, 
which only depens upon i/j and is the same for all nuclei: we then say that the 
RFG model superscales, namely / does not explicitly depend upon the momen- 
tum transfer q (scaling of the first kind) nor the Fermi momentum kp (scaling of 
the second kind). 

In the A-resonance peak, centered in uj = y/q 2 + mi — m^, all the above 
expressions still hold providing an appropriate scaling variable 




(fVp+Vr - - tun) , (7) 



accounting for inelasticity through the parameter p = 1 + (rn^ — m 2 N ) / (Armj^ ) , 
is used. The RFG model superscales also in the A region and the associated 
scaling function is the same as in the quasielastic peak. 
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3 The SuperScaling Analysis (SuSA) of electron scattering 

An extensive analysis of the world electron scattering data based on the su- 
perscaling ideas introduced above has been performed in the quasielastic do- 
main [10] by studying a reduced cross section, obtained as the ratio of the ex- 
perimental (e, e') cross section and the appropriate single-nucleon factors, as a 
function of the scaling variable (|6j. 

The analysis has shown that for high enough energies scaling of the first kind 
is fulfilled at the left of the quasielastic peak (the so-called scaling region) and 
broken at its right, whereas scaling of the second kind is very well satisfied at 
the left of the peak and not so badly violated at its right. Moreover, the scaling 
violations have been found to reside mainly in the transverse channel, whereas 
in the longitudinal scaling is very good. 

As a result a phenomenological superscaling function f(i/j) has been con- 
structed which fits the electron scattering longitudinal response for all the avail- 
able momentum transfers and nuclear species. The deviations from the total 
(longitudinal and transverse) data can be ascribed to two contributions: 1) the 
excitation of the A-resonance and 2) q- and k f -dependent nuclear effects such 
as meson exchange currents and their associated correlations, which break both 
kinds of scaling. Although the exact calculation of these effects is very in- 
volved [11-14], their net effect in the quasielastic region can be assessed to 
be roughly 10-15%. 

The superscaling function, plotted in Fig. 1 versus ip' (the "prime" indicates 
an energy shift, of about 20 MeV, necessary to reproduce the right experimental 
position of the quasielastic peak), significantly differs from the RFG one due 
to nuclear interactions not included in the Fermi gas model. In particluar, it 
is lower, it extends over a wider range of ip' and, importantly, it displays an 
asymmetric shape, with a pronounced tail in the i/>' > region. This peculiar 
feature of the superscaling function, which has been further investigated in re- 
cent work [2,5, 15], represents a stringent constraint on any reliable microscopic 
model. 

The same analysis has been recently extended to the A peak [8] by first sub- 
tracting the quasielastic cross section calculated by means of the phenomenolog- 
ical superscaling function from the total experimental data and then analyzing 
the result in terms of the scaling variable i/ja defined in Eq. @. This procedure 
removes from the data the impulsive (i.e., elastic eN) contributions that arise 
from quasielastic scattering. The corresponding reduced cross section (namely, 
divided by the appropriate elementary TV — > A function) has been shown to su- 
perscale reasonably well for negative values of ip&, indicating that the dominant 
process contributing to the cross section in this region has been correctly iden- 
tified as the N — ► A excitation. For positive values of tpA this is no longer the 
dominant process since other resonances and the tail of deep inelastic scattering 
come into play. 

As a result a phenomenological superscaling function / a (V>a)> fitting the 
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Figure 1. The experimental quasielastic superscaling function together with its phe- 
nomenological fit (dashed green) and the RFG result (solid red) plotted versus the corre- 
sponding scaling variable %p'. 

subtracted data, has been constructed in analogy with what was done in the 
quasielastic case (see Fig. 2). 
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Figure 2. The experimental A superscaling function for the world data on carbon and 
oxygen and its phenomenological fit (solid line) plotted versus the corresponding scaling 
variable ip' A . 

A test of the two phenomenological superscaling functions is presented in 
Fig. 3, where the data are compared with the result obtained by inserting these 
functions into the RFG responses 10} in place of /rfg- In Ref. [8] the same 
test is performed at different kinematical conditions and for different nuclei, 
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Figure 3. The experimental (e, e') cross section for 12 C compared with the SuSA predic- 
tion including both the quasielastic and A contribution (solid curve). The dashed curve 
represents the pure quasielastic cross section). 



4 Application of SuSA to neutrino and antineutrino scattering 

Having succeded in representing the inclusive electron scattering data by means 
of two "universal" (namely, q- and A-independent) superscaling functions, valid 
in the quasielastic and A-dominance regions respectively, we can now reverse 
the procedure and give predictions for neutrino and antineutrino reactions. 

Let us first focus on the charged-current process. Since the kinematics of 
(e, e') and (v, fi) are strictly related, the procedure simply amounts to mul- 
tiplying the phenomenological scaling function by the appropriate elementary 
neutrino-nucleon factor in order to construct the differential cross section (0. 

The assumption underlying this procedure is that the susperscaling function 
associated to the three responses in Q is the same and is equal to the one associ- 
ated to the electromagnetic responses in Q. While this is guaranteed by CVC as 
far as the vector-isovector components of the responses are concerned (i.e. the 
ones arising from the vector nuclear current), for the axial components it is not 
exactly true, since the axial current is only partially conserved. However, at the 
intermediate-to-high energies we are considering here, both the axial and vector 
operators reduce, in the non-relativistic limit, to er 1 : hence the corresponding 

1 Although we deal with relativistic effects exactly, this argument can be made more clearly in 
the non-relativistic limit. 
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scaling functions should be the same. This argument neglects higher order con- 
tributions like meson exchange currents, which of course act differently in the 
different channels, and go in any case beyond the present analysis. 

In Fig. 4 we show the double differential quasielastic CC neutrino cross sec- 
tion with respect to the outgoing muon momentum k' and solid angle O as a 
function of k' , for neutrino energy of 1 GeV and muon scattering angle of 90 
degrees. 

It appears that the superscaling analysis prediction (solid red) is significantly 
lower and wider that the RFG result (dashed green), reflecting the differences 
observed in Fig. 1. For comparison the non-relativistic Fermi gas result (dotted 
violet) is also shown in order to stress the importance of relativistic effects in 
this energy domain. Indeed these not only move the peak position and produce a 
shrinking of the response region due to the relativistic kinematics, but they also 
yield a change in size due to the correct Dirac "spinology". It clearly appears 
from Fig. 4 that relativistic effects cannot be neglected. Whereas they can be ac- 
counted for exactly in the simple Fermi gas model, some efficient prescriptions 
which can be used in order to "relativize" more sophisticated nuclear calcula- 
tions have been recently developed and successfully tested on the relativistic 
shell model [1]. 




100 150 200 250 300 350 400 450 500 550 600 
k'[MeV] 

Figure 4. Double differential (z/ M , fiT) cross section versus the outgoing muon momen- 
tum k' for neutrino energy E v = 1 GeV and muon scattering angle # M = 90° . Dashed 
(green): relativistic Fermi gas; dotted (violet): non-relativistic Fermi gas; solid (red): 
SuSA result. 

In Fig. 5 the neutrino and antineutrino differential cross sections are shown 
for E v — 1 Gev and scattering angle of 45 degrees, including also the A con- 
tribution (left peak). Note that the antineutrino cross section is about 5 times 
smaller than the neutrino one, the difference becoming more dramatic as the 
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scattering angle increases (see Ref. [8]). The reason is that the contributions 
of the responses Rt and R' T to the cross section are almost equal in size, but 
they add for neutrino and interfere destructively for antineutrino scattering. As a 
consequence small changes in the model, as the inclusion of two-body currents, 
could have very large effects on the 77 result, which should therefore be taken 
with great caution. 




300 400 600 800 1000 200 400 600 600 1000 



k'(MeV/c) k'(MeV/e) 

Figure 5. The SuSA cross section for the charged-current reactions (y^ , ^ ) (left panel) 
and (77 M , (right panel) on 12 C plotted versus the final-state muon momentum k'. The 
dash-dotted curves give the QE contribution, the dashed curves the A contribution and 
the solid curves the total. 

In the case of neutral current reactions, the kinematics are different from the 
one of CC processes and of (e, e'). In fact, while the CC reaction is a i-scattering 
type process (the Mandelstam variable t is fixed), the NC is a u-scattering pro- 
cess, since the detected final state is the outgoing nucleon, and the neutrino 
kinematic variables are integrated over [16]. This in turn implies an integration 
region in the residual nucleus variables which is different in the two cases. As 
a consequence, it is not obvious that the superscaling procedure, based on the 
analogy with inclusive electron scattering, is still valid. 

This issue is discussed at length in Ref. [9], where it is shown that the scaling 
method is based on a factorization assumption which has been tested numeri- 
cally, with the outcome that the procedure can be applied also to neutral current 
reactions. 

In Fig. 6 (left panel) we show the results for quasielastic neutrino scattering 
and proton knockout calculated in both the RFG model (dashed green) and the 
SuSA approach (solid red). More results, corresponding to antineutrino scatter- 
ing and to neutron knockout, are given in Ref. [9]. 

As for the CC case, it appears that the results based on superscaling are sig- 
nificantly different from the ones of the non-interacting model and the almost 
parabolic RFG cross section is strongly deformed by the nuclear dynamics im- 
plicit in the phenomenological scaling function. 

Finally, in Fig.6 (right panel) we show results obtained by allowing for a non- 
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zero strange quark content of the nucleon: we compare our predictions for the 
cross section in a situation where no strangeness is assumed (solid red) with the 
ones obtained including strangeness in the magnetic (dashed green) and axial- 
vector (dotted violet) form factors, using for fi s — G^ (0) a representative value 
extracted from the recent world studies of parity-violating electron scattering 
and taking g s A = G^\o) to be —0.2 [17]. The effects from inclusion of electric 

(s) 

strangeness are not shown here, since G E has almost no influence on the full 
cross sections. 

These results, which show a strong sensitivity of the cross section to the 
strangeness content of the nucleon, are in line with the well-known fact (see, 
e.g., Refs. [16, 18]) that the NC reactions can be used, together with parity- 
violating electron scattering [19,20], to measure the strange form factors of the 
nucleon. 




T r [MeV] T p \MeV] 



Figure 6. Double differential neutral current neutrino cross section for proton knockout 
with respect to the proton energy and solid angle plotted versus the proton kinetic energy. 
The neutrino energy energy is 1 GeV and the proton scattering angle 20 degrees. Left 
panel: the SuSA prediction (solid red) is compared to the RFG result (dashed green). 
Right panel: effects of variations of the magnetic (dashed green), axial (dotted violet) or 
both magnetic and axial (dotdashed cyan) strange form factors of the proton, compared 
with the result without strangeness (solid red). 



5 Conclusions 

We have shown how the large amount of existing (e, e') data can be used to 
predict neutrino-nucleus cross sections, of interest for the search of neutrino 
oscillations and for the measurements of the nucleon strange form factors. 

The method we have developed is based on the superscaling properties dis- 
played by the electron scattering data both in the quasielastic peak and in the 
A-excitation region. 

Beyond the practical usefulness of the SuSA procedure, this study opens 
several interesting questions connected to the microscopical description of the 
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superscaling function and in particular to its strongly asymmetric shape, which 
represents a powerful test of different nuclear models. 

We believe that, for energies not too low, the superscaling analysis yields 
predictions for neutrino scattering at the 15-20% level, the uncertainty being 
related to the violations of superscaling. These mainly arise from two body con- 
tributions (meson-exchange currents and their associated correlations) which are 
not included in our phenomenological representation of the nuclear dynamics. 
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